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We analyze and scattering with the microscopic coupled-channels method and investigate 

the coupled-channels and three-nucleon-force (3NF) effects on elastic and inelastic cross sections. In the mi¬ 
croscopic coupled-channels calculation, the Melbourne g-matrix interaction modified according to the chiral 
3NF effects is used. It is found that the coupled-channels and 3NF effects additively change both the elastic 
and inelastic cross sections. As a result, the coupled-channels calculation including the 3NF effects signifi¬ 
cantly improves the agreement between the theoretical results and the experimental data. The incident-energy 
dependence of the coupled-channels and 3NF effects is also discussed. 
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Introduction. Microscopic description of many-body nu¬ 
clear reactions is a long standing subject in nuclear physics. 
Elastic scattering, which is one of the most basic processes, 
was first understood phenomenologically with the optical 
model. The model was founded by the multiple scattering 
theory from the microscopic point of view. Accord¬ 

ing to the theory, elastic scattering is described by multistep 
processes with a nucleon-nucleon effective interaction. So 
far many types of density-dependent complex interactions {g- 
matrix interactions) have been proposed as effective interac¬ 
tions based on the Brueckner theory. It is, however, still not 
easy to reproduce even the elastic scattering data microscopi¬ 
cally. 

Recently, the Melbourne group achieved a great success to 
describe proton scattering with no adjustable parameter. They 
applied a p-matrix interaction constructed with the Bonn po¬ 
tential 101 to the folding model calculation, and reproduced 
the measured elastic cross sections and analyzing powers 
in a wide incident-energy range a. In addition, the Mel¬ 
bourne interaction was also applied to neutron scattering la 
and ®He-proton scattering 10,Si- Furthermore, total reaction 
cross sections of nucleus-nucleus scattering were well repro¬ 
duced sm. These are the remarkable progress of micro¬ 
scopic reaction theory in recent years. In some cases, how¬ 
ever, nucleus-nucleus elastic cross sections cannot be repro¬ 
duced d. One of the reasons for this failure may be three- 
nucleon force (3NF) effects on nucleus-nucleus scattering; the 
Melbourne interaction does not include the effects explicitly. 

3NF effects on many-body nuclear reactions are a hot topic 
relating to the equation-of-state of nuclear matter. In the g- 
matrix folding model, the 3NF effects are represented through 
the density dependence of g-matrix. In Refs. d[II, the re¬ 
pulsive nature due to 3NF effects is simulated by changing the 
two-nucleon force (2NF) in the nuclear medium to reproduce 
the empirical saturation. The g-matrix interaction (CEG07) 
obtained by including the “3NF” effects gives a less attractive 
and less absorptive potential than the original one. The ef- 
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fects are rather small for proton scattering but the agreement 
between the theoretical results and the experimental values 
are improved for the spin observables at intermediate ener¬ 
gies 111^ . For nucleus-nucleus scattering, the measured elas¬ 
tic cross sections even at large scattering angles are well re¬ 
produced ifoil . In Ref. ifldl] . as a first attempt of explicitly 
applying the phenomenological 3NF to nucleon optical po¬ 
tential, the Urbana 3NF llT^ and the three-nucleon interac¬ 
tion dill were used to construct 5-matrix interactions. The 
5-matrices were applied to proton elastic scattering, which 
showed the 3NF effects are again small. 

In the previous works lUu th e roles of 3NFs based on 
the chiral effective field theory 111^ I20I] were investigated in 
nucleon-nucleus and nucleus-nucleus elastic scatte ring . One 
of the advantages of using the chiral interactions 111914^ is 
that the two-, three-, and many-body forces are treated system¬ 
atically so that the uncertainty of interactions is minimized. At 
present, the chiral interactions are most reliable for investigat¬ 
ing the 3NF effects although they cannot be applied to high- 
energy reactions beyond the cutoff energy. In Refs, 
the Melbourne interaction was modified according to the ef¬ 
fects of the chiral NNFO 3NF based on the knowl¬ 

edge that 5-matrices from the Bonn potential and the chiral 
2NF are resembling; see Refs. |[Tll [3 for details. The chi¬ 
ral 3NF effects make the real (imaginary) part of the folding 
potential less attractive (much absorptive). This provides us a 
new insight into the 3NF effects, which had not been acquired 
in the studies employing phenomenological 3NFs. As a re¬ 
sult, the chiral 3NF effects improve the agreement between the 
theoretical calculation and the experimental data for nucleus- 
nucleus elastic scattering. The recent study by Toyokawa et 
al. 1I23I] in which g-matrices are constructed directly using the 
chiral 2NFs and 3NFs reinforces previous findings. 

Although the 3NF effects improve the description of 
nucleus-nucleus elastic scattering, there remain some discrep¬ 
ancies. The theoretical cross sections at large angles overshoot 
the data for scattering, whereas the theoretical result 

is slightly off the data in the diffraction pattern of 
scattering. One of the reasons for this discrepancy will be 
the coupled-channels effects due to the projectile and target 
excitations. In fact, these effects have been discussed in the 
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preceding studies 1^ based on the microscopic coupled- 
channels method. In the analyses inelastic cross sections were 
mainly focused and elastic cross sections were used just to 
determine the real and/or imaginary parts of the optical poten¬ 
tials. On the other hand, our interest in this report is to clarify 
the coupled-channels effects on both the elastic and inelastic 
nucleus-nucleus scattering. 

For this purpose, we will perform microscopic coupled- 
channels calculations, without any adjustable parameters, of 
i60-^®0 and scattering using the Melbourne inter¬ 

action modified with the chiral 3NF effects, and clarify the 
coupled-channels and 3NF effects on scattering observables. 
These two effects are treated simultaneously within the mi¬ 


croscopic coupled-channels framework, and their incident- 
energy dependence is also investigated. 

Framework of coupled-channels calculation. Since the mi¬ 
croscopic coupled-channels method to describe elastic and in¬ 
elastic scattering is well established, we give just a brief re¬ 
view of it; see, e.g.. Refs. ll26l - l^ for details. We consider 
scattering of a projectile P off a target T, where P and T are 
identical bosons. The radial part x^l partial 

wave between P and T, where L is the P-T relative orbital an¬ 
gular momentum and J is the total angular momentum, is ob¬ 
tained by solving the following coupled-channels equations. 
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Here, s — R rp — rrp, p = rp — s/2, and t = r^p s/2. 
The channel number indicating the intrinsic spins of P (Ip) 
and T (/rp) is specified by 7, and the channel spin S defined 
by |/p - Irp I < S < /p -f /rp. The subscript 0 indicates 
the initial channel. E..^ is defined hy E..^ = E — with E 
the incident energy in the center-of-mass system and e.^ the 
sum of the excitation energies of P and T. M is defined as 
M = ApAt/(A-p -f At), where Ap(T) is the mass number 
of P (T). U^°^^{R) is the Coulomb potential between P and T. 
W {LSL' S'\J\) is the Racah coefficient and = [2(1 -f 

<57p7p4psp)]-'/^- 

The input of the microscopic coupled-channels calculation 
is the nucleon-nucleon effective interaction and 
and the transition densities p^Aj, (A = P and T). In this 
paper, we use the Melbourne q-matrix interaction modified 
with the chiral 3NF effects Olll [Tst] . Although a g matrix con¬ 


tains an effect of single-particle excitation in nuclear matter, 
collective excitation for a specific nucleus is not adequately 
included. Following preceding studies, therefore, we take 
into account these collective excited states using the coupled- 
channels framework. A possible double-counting of the cou¬ 
pling to non-elastic channels is expected to be negligible, be¬ 
cause the channels explicitly included are specifically collec¬ 
tive and have no analogue in nuclear matter. 

We take the frozen density approximation for evaluating p 
in the argument of the g matrix; p = Pp{r^) -f prp{r^) is 
used, where means the midpoint of the interacting two 
nucleons. For the coupling potentials, we take the average 
of the densities in the initial and final states, i.e., = 

+ ptu)/2 for A = P or T. 

Since we need not only the ground state density but also 
transition densities for excited states, we adopt microscopic 
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cluster models to obtain them. For we consider the 0^, 
2 + (4.44 MeV), 0+ (7.65 MeV), and 2+ (10.3 MeV) states, 
and we use the transition densities obtained by the resonat¬ 
ing group method (RGM) based on a 3a model ll^ . These 
densities are found to reproduce the elastic and inelastic form 
factors determined by electron scattering and thus highly re¬ 
liable. For 1^0, the 0^", 3^ (6.13 MeV), and 2^ (6.92 MeV) 
states are taken into account. We use the transition densities 
obtained by the orthogonality condition model (OCM) assum¬ 
ing an structure ll^ . However, the transition densities 

between the excited states are not prepared; we thus neglect 
the coupling between the excited states in the calculation of 
the scattering. It is noted that we multiply the 

transition density between the 0^ and 3^ (2^) states by 0.816 
(1.17), so as to reproduce the experimental values of i?(E3) 
and B(E2), i.e., 1490 ± 70 fm® and 39.3 ± 1.6 fm'‘ 01, re¬ 
spectively. A similar procedure was adopted in Refs. 12611^ 
using the same transition densities of In the reaction cal¬ 
culation, we do not use the excitation energies of and 
evaluated by the microscopic models but adopt the experimen¬ 
tal values of them. 

In the present calculation, we treat transitions only by the 
nuclear interaction, though the Coulomb interaction is in¬ 
cluded in the diagonal components of the coupling poten¬ 
tials. The symmetrization of the identical bosons as well as 
the mutual excitation is included properly. We adopt the rel¬ 
ativistic kinematics of P and T. As for the integration for cal¬ 
culating the coupling potentials in Eq. Q, we perform the 
Monte-Carlo integration with random numbers generated by 
the Mersenne-Twister method 10 . 

Results and discussions. Figure [T] shows the differ¬ 
ential cross sections for the scattering at (a) 

70MeV/nucleon and (b) 44MeV/nucleon as a function of the 
scattering angle 6 in the center-of-mass system. In each panel, 
three cross sections corresponding to the 0]*' (top), 3j" (mid¬ 
dle), and 2^ (bottom) states of the ejectile are shown; the other 
particle is in the ground state in the final state. For the 0^ 
state, i.e., the elastic scattering, the ratio to the Rutherford 
cross section is plotted. The solid (dashed) line is the result 
of the coupled-channels calculation with (without) the 3NF 
effects, and the dotted (dot-dashed) line corresponds to the re¬ 
sult of the one-step calculation with (without^e 3NF effects. 
The experimental data are taken from Refs. 11341 - 1^ . 

From Fig. 1(a), we see that the calculation including both 
the 3NF and coupled-channels effects well reproduces the 
elastic scattering data up to 0 = 15°. The 3NF and coupled- 
channels effects are found to be comparable at forward an¬ 
gles 9 < 12°, whereas the former is more important than the 
latter at larger angles. As shown in Ref. fill] , the 3NF ef¬ 
fects change mainly the interior part of the nucleus-nucleus 
potential, which is influential to the scattering at larger an¬ 
gles. On the other hand, it is well known that in general the 
coupled-channels effects generate the so-called dynamical po¬ 
larization potential in the nuclear surface region. The separa¬ 
tion of the regions sensitive to the 3NF and coupled-channels 
effects makes the two effects rather independent. In fact, one 
sees that the 3NF (coupled-channels) effects do not change 
essentially whether the coupled-channels (3NF) effects are in- 




FIG. 1; (Color online) Differential cross sections of elastic and 
inelastic scattering at (a) 70 MeV/nucleon and (b) 44 

MeV/nucleon. The solid (dashed) line corresponds to the result with 
(without) the 3NF effects, and the dotted (dot-dashed) line corre¬ 
sponds to the result of the one-step calculation with (without) the 
3NF effects. The experimental data are taken from Refs, 

The inelastic cross sections are scaled by the factor shown inside the 
panel for clarity. 


eluded or not. 

For the 3^ and 2f inelastic cross sections, the solid line re¬ 
produces the data for 0 < 5° but slight overshooting is seen 
at larger angles. One sees that the 3NF effects improve the 
agreement with the inelastic cross section data. The 3NF ef¬ 
fects on the inelastic cross sections appear at large angles as 
in the elastic cross sections, although the effects are relatively 
small. In general, inelastic cross sections at very forward an¬ 
gles are mainly determined by the strength of the coupling 
potential. The negligibly small difference between the solid 
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and dashed lines at forward angles suggests the small 3NF 
effects on the coupling potentials. It is found that the 3NF ef¬ 
fects on the inelastic cross sections come from the changes in 
the diagonal part of the potential for each channel. It 

was pointed out in Ref. ll^ that the inelastic cross sections 
may be varied by couplings with even higher excited states. 
This can be a reason for the discrepancy for the 3j” and 2^ 
cross sections at larger angles and should be investigated in 
the future. 

At 44 MeV/nucleon (Fig. 1(b)), the coupled-channels ef¬ 
fects are larger and the 3NF effects are smaller than at 70 
MeV/nucleon. This is because that at lower energies the re¬ 
action takes place in the peripheral region. Nevertheless, the 
3NF effects are still non-negligible at large scattering angles. 
To be precise, the correspondence between the calculation and 
the data is worsened from at 70 MeV/nucleon. One of the rea¬ 
sons for this will be the fact that at lower incident energies the 
higher excited states, which are not taken into account, may 
play a more significant role. 

Figure 2 shows the differential cross sections for 
the scattering at (a) 85 MeV/nucleon and (b) 

30 MeV/nucleon. The meaning of the lines is the same as in 
Fig.m We show only the elastic and the 2^ inelastic cross sec¬ 
tion, for which the experimental data are available. It should 
be noted that the 0^, 2^, 0^, and 2^ states of are taken 
into account in the calculation as mentioned before. The ex¬ 
perimental data are taken from Refs. jjTl - f^ . 

At 85 MeV/nucleon, we find larger coupled-channels ef¬ 
fects than those in the scattering. As a result, the 

solid line decreases significantly and undershoots the data. It 
is found that the coupling between the 0^ and 2f states plays 
a dominant role in the scattering. Unfortunately, 

coupled-channels effects give no change in the diffraction pat¬ 
tern. For the inelastic cross section the solid line well agrees 
with the data. At 30 MeV/nucleon, as shown in Fig. 2(b), the 
3NF effects do not affect to both the elastic and inelastic cross 
sections. The measured elastic and inelastic cross sections are 
well reproduced by the coupled-channels calculations, except 
for the small difference at large scattering angles. 

In Table |D we show the total reaction cross sections ctr 
obtained by the coupled-channels and one-step calculations; 
the 3NF effects are included in both calculations. One sees 
that (7 r increases by 2-A% depending on the strength of the 
coupled-channels effects. For the scattering, the cal¬ 

culated results are slightly larger than the experimental val¬ 
ues iliflidd]] . This may be due to the fact that the root-mean- 
square (RMS) radius of the proton density calculated by the 
ROM is larger than its empirical value by 3 %. Since ctr is 
very sensitive to the proton and neutron RMS radii, some fine 
tuning of the microscopic density will be needed for a detailed 
analysis of ctr. 

Summary. We have analyzed the elastic and the 3^ and 
2^ inelastic cross sections for the scattering at 70 

and 44 MeV/nucleon, and also the elastic and the 2f in¬ 
elastic cross section for the scattering at 85 and 30 

MeV/nucleon, by means of the microscopic coupled-channels 
method. The coupled-channels and 3NF effects on these cross 
sections were investigated. 



FIG. 2: (Color online) Differential cross sections of elastic and 
inelastic scattering at (a) 85 MeV/nucleon and (b) 30 

MeV/nucleon. The meaning of the lines is the same as in Fig. [T] 
The experimental data are taken from Refs. 


Since the coupled-channels and 3NF effects additively 
change both the elastic and inelastic cross sections, it is im¬ 
portant to consider the two effects together for a quantitative 
comparison between the calculation and the experimental data 
for nucleus-nucleus scattering. It was found that the 3NF ef¬ 
fects on the inelastic cross sections are almost the same as on 
the elastic cross sections. 

In detail, the relative importance of the 3NF and coupled- 
channels effects depends on the incident energies and scatter¬ 
ing angles. A key feature for understanding this dependence is 
the fact that the 3NF effects strongly affect the interior part of 
the nucleus-nucleus potential, whereas the coupled-channels 
effects affect the surface part of it. As a general tendency, we 
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TABLE I: Calculated and measured total reaction cross sections ctr 
( in the unit of mb). The results calculated with and without coupled- 
channels effects are shown; the 3NF effects are taken into account in 
both calculations. 


system 

CC 

Ich 

Exp. 

O- O @ 70Me V/nucleon 

1453 

1426 

- 

0- 0 @44Me V/nucleon 

1519 

1473 

- 

i2C-^^C@85Me V/nucleon 

1078 

1052 

998 ± 13“, 960 ± 30* 

i=C-^^C@30Me V/nucleon 

1262 

1214 

1209 ± 32“, 1315 ± 40* 


“At 86.3 MeV/nucleon na 
*’At 83 MeV/nucleon ITlll 
“At 32.5 MeV/nucleon 1401 
^At 30 MeV/nucleon S 

can say that the 3NF (coupled-channels) effects are important 
at relatively high (low) incident energies and at large (small) 
scattering angles. It was found that the 3NF effects are negli¬ 
gibly small at lower than about 30 MeV/nucleon. 

We have also shown that the 3NF effects on the coupling 
potentials are quite small. This may allow one to neglect the 


3NF effects in the analysis of the inelastic cross sections at 
very forward angles. The coupled-channels effects are found 
to increase the total reaction cross sections by about a few 
percent. 

In conclusion, the coupled-channels calculation including 
the 3NF effects significantly improves the agreement between 
the theoretical results and the experimental data. However, 
there remain some discrepancies between them, at backward 
angles in particular. It will be important to perform coupled- 
channels calculations including even higher excited states. A 
more detailed examination of the effective interactions and 
transition densities will be another important subject. 
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